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1  Introduction 


Two  of  the  profound  open  problems  in  the  theory  of  three  dimensional  viscous  flow  are 
the  unique  solvability  theorem  for  all  time  and  the  existence  theorem  for  the  global  attrac¬ 
tor.  We  have  shown  in  our  earlier  studies  [10,  9,  11,  13]  that  certain  regularizations  of  the 
Navier-Stokes  equations  are  uniquely  solvable  (up  to  dimension  six)  and  can  be  character¬ 
ized  by  compact  global  attractors.  A  natural  question  then  is  to  investigate  the  possibility 
of  establishing  such  results  for  the  conventional  Navier-Stokes  equations  by  a  limit  process. 

In  this  paper  we  will  accomplish  this  for  the  two  dimensional  case.  We  prove  in  particular 
that  the  compact  global  attractor  Ae  for  the  regularized  system  converges  to  the  compact 
global  attractor  A  of  the  conventional  system  as  e  — >  0. 

The  outline  of  the  paper  is  as  follows.  In  §2  we  present  the  relevant  mathematical  frame¬ 
work  for  the  paper.  In  §3  we  establish  several  uniform  estimates  for  the  generalized  solutions 
of  the  regularized  system.  These  estimates  hold  independent  of  the  size  of  the  regulariza¬ 
tion  parameter  and  remain  valid  when  this  parameter  goes  to  zero.  In  §4  we  show  that  the 
solution  of  the  regularized  system  converges  to  the  solution  of  the  conventional  system  as 
the  regularization  parameter  goes  to  zero.  The  next  two  sections  deal  with  convergence  of 
attractors.  In  §5  we  establish  the  convergence  of  the  simple  case  of  time  periodic  solutions. 

Then  in  §6  we  prove  the  central  result  of  this  paper  establishing  the  upper  semicontinuity  of 
the  compact  global  attractor  Ae  at  e  =  0. 

The  question  of  the  convergence  of  Ae  to  A  is  thus  completely  answered  for  the  two  di¬ 
mensional  case.  Our  future  investigations  will  be  concerned  with  the  corresponding  problem 

for  the  three  dimensional  case.  We  have  already  reported  partial  results  in  this  direction  in - - 

r 

[10,  9]  where  it  was  shown  that  the  generalized  solutions  of  the  three  dimensional  regular¬ 
ized  system  converges  to  the  solution  of  the  conventional  system  under  certain  conditions. 

:i 

The  methods  validated  in  this  paper  certainly  give  us  guidelines  to  elaborate  the  three  di-  - 

mensional  case.  We  remark  here  that  the  use  of  the  artificial  viscosity  method  to  establish  ^ 

solvability  has  also  been  successful  in  other  branches  of  partial  differential  equations.  A  well  •  Codes 

,  uid/or 


IHi.l  1  al 


□  □ 


known  example  of  such  a  result  is  the  viscosity  solution  method  for  the  Hamilton-Jacobi 
equations  [2]. 

2  Notations  and  Preliminaries 

We  regularize  the  Navier-Stokes  system  by  adding  a  fourth  order  artificial  viscosity  term 
(Laplacian  square)  to  the  conventional  system.  In  this  paper  we  will  restrict  ourselves  to 
periodic  boundary  conditions.  A  thorough  study  of  the  regularized  system  with  this  and 
other  types  of  boundary  conditions  was  carried  out  in  [10,  9,  11]. 

Let  us  consider  the  velocity  field  u  =  (ui,u2)  and  pressure  field  p  with  space  periodic 
condition  in  R2  such  that  for  0  =  (0,  L)  x  (0,  L),  we  have 


— — h  eA2u  —  i/A u  +  (u 
at 

•  V)u  +  Vp  =  f , 

in 

fi  x  (0,  oo), 

(1) 

<1 

e 

II 

o 

in 

U  x  (0,  oo), 

(2) 

u(x  +  Lei,t)  =  u(x,i) 

i  =  l,2, 

V 

t  6  (0,  oo), 

(3) 

p(x  +  Le{,t)  =  p(x,t) 

i  =  1,2, 

V 

t  e  (0,  oo), 

(4) 

u(*,0)  =  u0(*), 

in 

ft. 

(5) 

Here  v  >  0  is  the  coefficient  of  the  kinematic  viscosity  of  the  fluid  and  e  >  0  is  the  artificial 
dissipation  parameter.  /  and  uo  are  prescribed  vectorfields.  Notice  that  the  conventional 
Navier-Stokes  equations  correspond  to  the  value  of  the  regularized  parameter  e  =  0. 

We  denote  by  £Tm(fl),  the  Sobolev  space  of  L-periodic  generalized  functions  (condition 
(3))  which  have  up  to  order  m  square  integrable  distributional  derivatives.  These  spaces  are 
endowed  with  the  inner  product 

(u,  v)m  =  (Dau,Da  v)l»(o) 

|oc  |<m 

and  the  norm 

llullm  =  (  52  ll-c>aulli,*(fi))1/2- 

|Ot  |<m 
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Suppose  we  expand  «  G  Hm(£l)  by  the  Fourier  series 

«(*;  =  c*e  > 
kez* 

then  the  Sobolev  space  u  €  Hm(Q )  can  be  characterized  by 

=  {„|  c*  =  c_*,  £  l*l2mkl2  <  °°}- 

kez 3 

We  denote  by  Hm(Ct)  the  subspace  of  Hm(Q)  with  zero  average: 

Hm( fi)  =  {ti  €  ;  /  u(*)d*  =  0}. 

For  m  =  0,  we  have  H°(Q)  =  i,2(f2). 

We  introduce  the  following  solenoidal  subspaces  which  are  important  to  our  analysis: 

H  =  { V  e  L2(fi),  divu  =0,  U  •  n|r{  =  -U  •  n|ri+a,  i  =  1,2}; 

V  =  {u  G  H1^),  divu  =  0,  7ou|iv  =  7o«|ri+a,  *  =  1,2}. 

V  =  {«  €  i?2(n),  divu  =  0,  7o« |r<  =  7ou|ri+a,7i«lri  =  -7i«|r<+»,  *=  1,2}. 

Here  the  faces  of  fi  are  numbered  as 

Ti  =  dQ  fl  =  0}  and  I\+2  =  9f2  fl  {x*  =  L},  i  =  1,2. 

Here  7o,7i  are  the  trace  operators  and  n  is  the  unit  outward  normal  on  dfl. 

The  space  H  is  endowed  with  the  inner  product  (u,r)£,»(n)  and  norm  |u|  =  lu>w)I*(n) 
(or  ||  •  ||o).  It  can  be  shown  that  the  norm  induced  by  Hx{Tl)  and  the  norm  ||Vu||£,a(n) 
are  equivalent  in  V.  Similarly,  in  V  the  norm  induced  by  H2(Tt)  is  equivalent  to  the  norm 
||AuIU*(n)-  We  then  denote  ||u||  =  ||Au||£,a(n)  =  (u,u)y2  as  norm  m  V  derived  from  the 
inner  product 

(u,v)y  =  (Au,  At>). 

Let  V',  V'  denote  the  dual  spaces  of  V  and  V  respectively. 
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The  Stokes  operator  Ai  can  be  characterized  explicitly  using  Fourier  series.  We  write  the 


Stokes  problem, 


—  Au  +  V?  —  g,  in  Q, 


V-t*  =0, 


in  fi. 


We  have  in  fact, 


4i«  =  g,  with  D(Ai)  =  {u  G  H,  A\u  G  H}. 


«  G  D(Ai)  =  H2{Q)  n  H  =  V. 


Similarly,  we  can  solve  by  Fourier  series  the  following  linear  problem  which  is  fundamental 


to  our  analysis: 


A2v  +  Vp  =  /  ,  in  n, 


[  V  ■  v  =  0,  in  f2. 

The  linear  operator  A  (which  we  call  the  dissipation  operator)  is  characterized  by  Av  =  / .  It 
can  be  shown  that  the  self  adjoint  operator  A  G  C(D(A)\  H)  D  C(V;V')  is  closed  with  D(A) 
dense  in  V  C  H.  Moreover,  we  can  define  positive  as  well  as  negative  powers  A“,a  G  R 
with  domain  D(Aa).  Let  us  denote  by  Xa  the  space, 


Xa  =  D(Aa/A)  =  {»  G  Ha( fi),  div«  =  0}. 

In  fact  the  norm  induced  by  Ha(Q)  is  equivalent  to  the  norm  |A“^4u|  in  Xa.  This  means 
AIMk<|^a«|</?2|H|4a,  VueD(Aa),  VaGfi.  ({ 

We  can  deduce  from  a  theorem  due  to  Lions  [7]  that  D{A1^2)  =  V.  Thus, 

D{Ai)  =  V  =  Z?(A1/2).  0 


By  Rellichr8  Lemma  [1]  A  1  as  a  mapping  in  H  is  compact.  Hence  the  spectrum  of 
operator  A  is  discrete  with  finite  multiplicities  and  can  be  written  explicitly  using  Fourier 


series  as 


Ifiir4 


4 


The  self  adjoint  operator  A  possesses  an  orthonormal  set  of  eigenfunctions  {tu  complete 
in  H, 


Awj  =  Xjio  j  ,  t i>j  £  D(A ),  V  j. 

The  above  results  can  be  applied  to  the  Stokes  operator  A\  as  well.  The  eigenvalues  of  At 
are 

The  corresponding  orthonormal  set  of  eigenfunctions  {ti  j}”  x  is  complete  in  H  and 

Axwj  =  A jWj  ,  wj  E  D(AX),  V j. 

Remark:  For  the  space  periodic  case  considered  in  this  paper,  the  Stokes  operator  At  is 
in  fact  equal  to  A1/3.  However  since  for  other  types  of  boundary  conditions  they  are  different 
we  prefer  to  give  them  distinct  notations. 

Let  us  now  define  the  trilinear  form  &(-,-,•)  associated  with  the  inertia  terms: 

b(u,v,w)  =  ^2  /  uiDiVjWjdx ,  Di~—. 

i,r=  i Ja  dx< 

It  can  be  easily  shown  that 

b(u,v,u>)  =  V  u,v,w  €  V 

and 

b(v,v,  w)  =  0,  V  u,v  G  V. 

In  the  case  of  periodic  boundary  conditions,  we  have  in  addition  [15,  Lemma  3.1] 

&(«,«,  At  «)  =  0,  VuGP(At).  (10) 

Using  the  discrete  Holder  inequality  and  Sobolev  embedding  theorem  we  can  show  that 
fe(- ,  -,  -)  is  trilinear  continuous  on  x  Hm,+1(Q )  x  >  0  [15]: 

|6(tt,»,w)|  <  Co||tt||milhl|m,+t||w|lm„ 

5 

i 


V 


mi  +  m2  +  m3  >  1. 


The  following  well  known  estimate  will  be  used  later, 

|  &(tx ,  u ,  to )  |  <  c1|«|1/2|Vttj1/2|t,|1/2|Vt,|1/2|Vw|,  G  V{  or  V ),  (11) 

where  Ci  is  a  positive  constant. 

The  continuity  property  of  the  trilinear  form  enables  us  to  define  (using  Riesz  represen¬ 
tation  theorem)  a  bilinear  continuous  operator  B  from  Hmi (fi)  x  Hmt+1(Cl)  into  (i/ms(Q))'. 
In  particular,  for  u,v,w  G  V,  B(u,v)  G  V'  will  be  defined  by 

<  J3(u,w),t0  >v'xv=  &(«,«, to),  V  to  G  V. 

Similarly,  we  define  f?(u,u  )  €  V'  by 


<  B(u,v),w  >g,x$=  b( u,v,w), 


Vto  G  V. 


Using  the  operators  defined  above  we  can  write  the  regularized  system  (l)-(5)  in  the 
evolution  form: 


+  eAv€  + vA\ut  + B(uttut)  =  /,  t  >  0, 

ue(0)  =  urf. 


(12) 


The  existence  and  uniqueness  theorems  for  initial  value  problem  (12)  can  be  found  in  [10,  9]. 
The  main  result  in  this  work  is:  . 


Theorem  2.1 

(i)  Let  f  G  0,  T\  V ')  and  tio  G  H  be  given.  Then  there  exists  a  unique  weak  solution  of 
(12)  which  satisfies  ut  €  C([0,T];.ff)  H  2/2(0, T\  V),  VT  >  0. 

(ii)  Let  f  G  L3(Q,T]H)  and  Uo  EV  be  given.  Then  there  exists  a  unique  strong  solution  of 
(12)  which  satisfies  ue  G  C7([0,  T}\  V)  D  £2(0,  T\  D(A))}  VT  >  0. 

□ 
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Let  us  denote  by  We(t)  the  nonlinear  semigroup  associated  with  the  solution  of  the  regularized 
system  (12).  This  means  We(t,  Ojueo)  =  uc(t)  for  all  t  >  0. 

Notice  that  the  conventional  Navier-Stokes  system  can  be  written  in  the  evolution  form: 


CL  U  *  ,  v 

—  -|-  vA\U  +  B(u,u)  =  /,  t  >  0, 
at 

u(0)  =  UQ. 


(13) 


The  following  unique  solvability  theorem  for  the  system  (13)  is  well  known. 


Theorem  2.2 

(i)  Let  f  £  L2(0,T]V')  and  u0  £  H  be  given.  Then  there  exists  a  unique  weak  solution  of 
(IS)  which  satisfies  u  £  C([0,T]‘,  H)  D  £2(0,T;  F),  VT  >  0  (Lions  and  Prodi  [8]). 

(ii)  Let  f  £  L2(0,T;  H )  and  u0  £  V  be  given.  Then  there  exists  a  unique  strong  solution  of 
(IS)  which  satisfies  u  £  CQO,  T)\  V)  D  L2(0,  T;  D(v4i)),  VT  >  0  (Kiselev  and  Ladyzhenskaya 

[5])- 


□ 


We  will  denote  by  W{t)  :  u0  — ►  «(t)  the  nonlinear  semigroup  associated  with  the  solution 
of  Navier-Stokes  equations. 

Let  us  introduce  the  following  result  [10,  9]  which  is  needed  for  the  convergence  of  ut  as 
the  regularized  parameter  e-+0. 

Proposition  2.1  Let  {ujt}  be  a  sequence  of  vectorfields  that  converge  weakly  in  L2(0,  T;  V), 
weak-Star  in  L°°( 0,  T;  H )  and  strongly  in  L2( 0,  T;  H).  Then,  the  following  limits  are  obtained 
for  any  vector  function  w  in  y  =  {u?  £  C{0,T',V)]w'  £  L3(0,T]  H)} 

(a)  lim  f  (uk(t),w\t))dt  =  f  (u (t),w\t))dt, 

fc— » oo  Jo  Jo 

(b)  lim  /T(Vu k(t),Vw(t))dt=  [T(Vu(t),Vw(t))dt, 

k—*oo  JO  Jo 

(c)  lim  /  b(uk(t),uk(t),w(t))dt=  I  b(u(t),u(t),w(t))dt. 

fc— >oo  Jo  Jo 

□ 
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3  Uniform  Estimates  for  the  Solutions 

In  this  section  we  will  establish  various  estimates  uniform  in  e  for  the  solution  of  reg¬ 
ularized  Navier-Stokes  equations.  These  bounds  will  be  used  to  establish  the  limit  of  these 
solutions  to  the  conventional  Navier-Stokes  equations. 

Lemma  3.1  Let  f  £  L2(0,T;  if),  then  the  weak  solution  u e(i)  of  the  regularized  Navier- 
Stokes  equations  satisfy 


(i)  sup  |u£(i)|2  <  du 

t€[0,T] 

(ii)  j* \Vue{t)\2dt  <  d2, 

where  di,^  are  both  independent  of  e. 


Proof:  By  taking  the  inner  product  of  (12)  with  ue(t),  we  obtain  the  energy  equality: 

+  2fHu.ll1  +  2r|Vu,|J  =  2 :(/,«,). 

Here  we  have  used  the  fact  that  b(ue,  ut,  ue)  =  0.  By  applying  Young’s  inequality  and  the 
Poincare  Lemma  |Vue|2  >  ni\ue\2,  we  get 

4l».|’  +  2f||..||1  +  u|V».|J  <  kJl,  '  (14) 

dt  u  ji  i 

where  fii  =  47r2 / L2  is  the  smallest  eigenvalue  of  the  Stokes  operator  A\.  If  we  drop  the 
positive  term  associated  with  e,  we  obtain 

<  kL. 

Hence,  by  integrating  the  above  inequality  from  0  to  t,  we  get 

MOP  < 

Ct  Jo 


with  a  =  vp.\.  That  is, 


sup  |ue(<)|2  <  d\. 

t€[0,T] 
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1  tT 

with  di  =  di(ut0t  f  ,a,T)  =  \ucol2  +  -  \f(s)\2ds. 

a  Jo 

If  we  drop  the  term  2e||«£||2  in  (14)  and  integrate  from  0  to  T , 

l»,(T)|!  +  V  /T|V««|Vf  <  +  i  [T\f\2dt. 

Jo  a  Jo 

This  implies 

fT\Vu.\2dt<-{  l-al'  +  i  [T\f\2dt], 

Jo  v  a.  Jo 

That  is 

/T|Vu£(i)|2di  <  d2,  d2  =  ^. 

Jo  v 

Notice  that  d\  and  d2  do  not  depend  on  e. 


(15) 


□ 


Corollary  3.1  Let  f  GE  H  be  independent  of  time  then, 

<  |»,(0)|V“‘  +  fil(  1  -  e-“).  Po  =  V'  v  ‘  >  °’ 

(ii)  fT  |V«.| 2dt  <  i(  |»«|J  +  !^),  T  >  0. 

Jo  v  a 

Proof:  By  integrating  (14)  from  0  to  t  with  /  G  we  obtain  result  (i)  with  p0  =  \f\/ot. 
Result  (ii)  is  a  direct  consequence  of  (15). 

□ 

Lemma  3.2  Let  f  €  L2(0,T\  H),  then  the  strong  solution  ut(t)  of  the  regularized  Navier- 
Stokes  equations  satisfy 

(i)  sup  |Vue(0|2  <  <*3, 

te[0,T] 

(ii)  iT\A,uc{t)\2dt<d<, 

Jo 

where  d3,dA  are  both  independent  of  e. 
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Proof:  We  consider  the  strong  solutions  and  take  the  inner  product  of  (12)  with  Ai ut 
obtaining 

jt\ V«e|2  +  2e|v43(/4ue|2  +  2u\A1ut\2  =  2  (/ ,  A1u€). 

Here  we  have  used  the  fact  that  b(ue,  ue,  Axut)  —  0  and 

(Au„  A,..)  =  A'**,) 

=  |  A'I'utf 

since  D(A\)  =  D(A1^2)  and  both  A,A\  are  self  adjoint.  By  applying  Young’s  inequality,  we 
get 

i Vue|2  +  2e |A3/4ue|2  +  v\Axut\2  <  1^.  (16) 

at  V 

We  can  drop  the  positive  term  2e  |v43/4ue|2  and  use  the  inequality  |>liue  |2>^|V  u€  | 2  to  get 

+  KfnlVu.l1  <  — .  (17) 

at  v 

Integrating  the  above  inequality  from  0  to  t,  we  get 

|V„,(i)|*  <  |V»<0pe""+  1  A/WI2 

v  Jo 

where  a  =  v\i\.  This  gives 

sup  |Vue(t)|2  <  d3. 

<€[0,T] 

with  d3  =  d3(ut0,f,u,T)  =  IVueol2  +  ~  /  |/(«)|3<£.s. 

V  Jo 

We  can  drop  the  term  2e  \A2^ut\2  in  (16)  and  integrate  from  0  to  T  to  get 

| Vut(T)|2  +  u  [T \A,ur\2dt  <  IVueol2  +  -  [T\f  \2dt.  (18) 

Jo  v  Jo 

This  implies 

F  | A^.I’Jt  <  i[  |V»«|!  +  I  fT  \f\2dt  ]. 

Jo  V  V  Jo 

That  is 

f  \Aiut(t)\2dt  <  (£4,  di  =  — . 

Jo  v 

Notice  that  d3  and  <£4  do  not  depend  on  e. 
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□ 


From  this  proof  we  can  easily  deduce  the  following  results. 

Corollary  3.2  Let  f  £  H  be  independent  of  time  then  the  strong  solution  ue(t)  satisfies: 

0)1  <  iv«<(o)i v»‘ + pK  i  -  *-“•),  =  kli  v  t  >  o, 

va 

rT  1  \f\2T 

(ii)  /  \Aiuc\2dt  <  — (  |Vueo|2  +  — — ),  T  >  0. 

Jo  V  V 

□ 

4  Limit  to  the  Navier-Stokes  Equations 

In  this  section  we  will  establish  the  convergence  of  the  strong  solution  of  the  regularized 
Navier-Stokes  equations  as  e  — »  0.  In  [10,  9],  a  similar  convergence  result  for  the  weak  solution 
was  established  under  certain  conditions.  In  this  section  we  will  prove  the  convergence 
without  any  assumption  on  the  bound  for  the  solutions. 

Theorem  4.1  Let  u €{t)  be  the  strong  solution  given  by  Theorem  2.1.  Then  as  e  — ►  0  , 
the  solution  u€  converges  to  a  unique  limit  which  is  the  strong  solution  of  the  Navier-Stokes 
equations. 

Proof:  We  need  three  forms  of  convergence  for  appropriate  subsequences.  Namely, 

(i)  u€h,  — >u  in  L2(0,  T;V)  weakly; 

(ii)  u tkl  — in  L°°(0,T-,H)  weak-star;  (19) 

(iii)  uefc,  — ►  t*  in  L2(0,T]H)  strongly, 

as  ek>  — +  0. 

It  follows  from  Lemma  3.2  that  utk  G  L2(0,T]  D(AX))  D  L°°(0,T;  V)  with  bounds  in¬ 
dependent  of  ek.  This  easily  implies  (i)  and  (ii).  (In  fact,  better  convergence  results  are 
obtained.)  The  strong  convergen''®  result  (iii)  can  be  established  as  follows: 
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Let  {«*,,}  and  {u£,}  be  two  sequences  of  strong  solutions  of  (12)  corresponding  to  the 
same  initial  data  with  e*  =  £  and  ej  =  respectively,  (i.e.  e*,  ej  — ►  0  as  k,  l  — ►  oo.)  We 
denote  by  t ie(t)  =  u£Ai(f)  —  ti£( (t)  so  that  u£(0)  =  0.  By  subtracting  the  equation  (12)  for  tt£( 
from  that  for  u£(t,  we  get 

u'  +  ekAuCk  -  etAuei  +  vAxue  +  B[uek,  u£J  -  B(u£l,  u£()  =  0  (20) 

Taking  the  inner  product  of  (20)  with  u£,  we  get 

(«e.«£)  +  <*(>!«<*,  U£)  -  tl(Auei,Ue)  +  v{A1Uc,U€ )  +  b{u£k,u€h,u£)  ~  6(«£, ,  U£, ,  U£)  =  0. 

By  using  the  fact  that  ( Au,v  )  =  {A,  u,Aiv)  and  the  properties  of  the  trilinear  form  b 
established  earlier,  we  obtain 

+  ejt(>li««*,^itie)  -  e/(Aiu£(,,4iu£)  +  i/|Vu£|2  -  6(u£,  «£,«£J  =  0. 

Thus 

|t*e|2  +  v\ Vm£|2  <  e*  |(^4iuefc,  i4itte)|  4-  ei  |(i4i«£|,  Ai«£)|  -f  |i>(u£,  «c,  «*»)!•  (21) 

The  trilinear  term  can  be  estimated  using  (11)  as, 

|fc(t*£,Ue,U£J|  <  C!  |ttc|  |Vu£|  |V«£J 

<^lV-eia  +  |;ittel2|V«eJ2. 

Substituting  the  above  result  into  (21),  we  obtain 

iu£|2  +  i/|Vu£|2  <  2e*  | A\ u €k J  \Axut\  +  2 ei  |j4iue|  +  —  |u£|2  |Vu£J2. 

at  v 

We  drop  the  positive  term  i/|Vu£|2  and  integrate  this  equation  from  0  to  t,  and  noting  that 

u  £(0)  =  0,  we  get 

|«£(i)|a  <  2e*  /  |i4iu£)k|  \A\ut\dr  +  2ej  [  |Aiut||  \Axu€\dT 
Jo  Jo 

(22) 

u  Jo 
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From  Lemma  3.2,  we  have  uflt  €  L2(0,T;  D(Ai ))  with  bound  independent  of  e*.  Hence  by 
applying  the  Schwartz  inequality  we  have 

£kfo  \Aiut\dr  <  ek(J^  |A1ueJ2dr)1/2(^  |i4lU e|2<£r)1/2 

<  £k(31/2d4). 


Similarly  we  have 


et  f  |Ai«ei|  \Axu€\dT  <  e/(3 1/2d4). 
Jo 


Since  from  Lemma  3.2  ueit  €  L°°(0,T ;  V)  with  bound  independent  of  e*, 


/  |««|2  |VueJ2dr  <  d3  I  |ue|2dT. 
Jo  Jo 


Combining  all  of  the  above  results  into  inequality  (22)  we  get, 


|u4(f)|2  <  2 . 3 1/Jd4(cfc  +  e,)  +  ^  r  |«e(r)| 2 dr.  (23) 

v  Jo 

Notice  that  here  C\,d3  and  d4  are  all  independent  of  e.  Now,  if  we  denote  y(i)  =  Jq  |ue  (r)|2dr 
then  (23)  becomes 

j  §•  <  2  ■  3 +  €,)  +  ^v(() 


y(0)  =  o. 


From  this  we  can  deduce  that, 


Hence, 


That  is 


y(i)  <  31/2i/(efc  +  e/)^  eXp  ~  1  ,  t6[0,T]. 

tT 

/  |ue(r)|2dr  — >  0,  as  efc,  6/  -►  0. 

Jo 


ueh  -  «ciI|l3(0,T;H)  - *  0,  as  €k,  6/  -*  0. 


This  means  {uefc}  is  a  Cauchy  sequence  in  L2(0,  T;  H)  and  there  exists  a  limit  u  in  L2(0,  T;  H) 


such  that 


“o,  —  «|UJ(o,Titf)  — *  0  as  efc  -+  0. 
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This  completes  the  proof  of  the  strong  convergence  result  (iii)  of  (19). 

Let  us  now  establish  the  limit  of  the  equation  (12)  as  ek  — >  0.  Let  us  take  a  test  element 
4>  such  that 

<f>  E  <7(0,  T\  V)  and  <j> '  e  L2(0 ,  T ;  H ). 

Taking  the  inner  product  with  (12),  we  get 


(u€„)^)  +  ek(AuCk,«f>)  +  i/(Aiuek,<f>)  +  b(uek,u€k)(l>)  =  (/,<£), 

Integrating  with  respect  to  t  and  then  integrating  by  parts  we  get, 

(u e„(n<t>(T))  -  (t^(0),*(0))  -  fT(uek,*')dt 

Jo 

rT  rT 

+£fc  /  (Aiutk,  Aitf>)dt  +  v  /  (Vu£fc,V0)dt 

Jo  Jo 

rT  rT 

+  Ku€k,^tk,d>)dt=  ( f,<t>)dt .  (24) 

JO  Jo 

Let  us  choose  V>  =  <t>  €  X>(0,T;  V).  Using  the  convergence  Proposition  2.1  we  can  take 
the  equation  (24)  to  the  limit  as  e*  — »  0, 

—  /  (« ,ij>')dt  +  i/  I  (Vu,  Vi/>)df -+•  /  b(u,u,rl>)dt  =  f  (f,ij>)dt. 

Jo  Jo  Jo  Jo 


Here  the  term 


as  ejt  — >  0  since 


rT 

*k  ( A1utk,Aiip)dt — ►  0, 

Jo 


J'  (Aiutk,  Air/}  )dt  <  const,  independent  of  e*. 
o 

Since  u  €  L2(0,T;jD(^))nL°°(0,T;  V),  we  can  conclude  that  u  is  indeed  the  strong  solution 
for  the  conventional  Navier-Stokes  equations. 


□ 

From  Lemma  3.1,  we  can  deduce  that  the  weak  solutions  uek  6  L2{ 0,  T;  V")  D  L°°(0,  T\  H) 
are  bounded  independent  of  e*.  Hence  the  convergence  results  (i)  and  (ii)  are  satisfied  by 
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the  weak  solutions.  However,  such  convergences  are  not  sufficient  to  conclude  the  strong 
convergence  in  L2(0,T;  H).  Therefore,  Theorem  4.1  is  not  valid  for  weak  solutions.  We  can 
nevertheless  establish  a  convergence  theorem  for  the  weak  solutions  provided  they  satisfy  a 
certain  bound.  The  following  theorem  is  proved  in  [10,  9]: 

Theorem  4.2  Let  u€  be  the  weak  solution  of  problem  (12)  given  by  Theorem  2.1  with  f  E 
L2(0,T]  H).  Moreover,  ue  E  I»°°(0,  T;  L3+S(Q))  uniformly  in  e  (6  >  0)  and  v  >  c,  where  c 
is  a  positive  constant.  Then ,  for  e  — >  0,  ue  approaches  a  unique  limit  u  which  is  the  weak 
solution  of  the  Navier-Stokes  equations. 

□ 

5  Convergence  of  T-time  Periodic  Solution 

In  this  section  we  will  show  that  the  T-time  periodic  solution  of  the  regularized  Navier- 
Stokes  equations  converges  in  the  limit  to  a  T-time  periodic  solution  of  Navier-Stokes  equa¬ 
tions.  Existence  theorems  for  T-time  periodic  solutions  to  the  regularized  Navier-Stokes 
equations  were  established  in  [10,  9].  For  similar  existence  theorems  for  conventional  Navier- 
Stokes  equations  see  [12,  17,  14].  Let  us  first  give  a  definition  for  the  time  periodic  solution. 

Definition  5.1  Let  f(t)  E  L2oe(-oo,-\-oo,,H')  and  be  periodic  with  period  T.  Then  «e  (0 
is  said  to  be  a  T - time  periodic  solution  of  regularized  Navier-Stokes  equations  if  it  satisfies 
ue(t)  G  L°°(—oo,  -foo;  H )  D  L2oc(  — oo,  -foo;  V),  T-periodic  in  time  and 

f  {-(««(<),  V»'(*))  +  e(Au*(t),  AV-(t))  +  KV«t(<),VV»(t)) 

J—oo 

+  b(u€(t),uc(t),rl>(t))}dt  =  f  (f(t),rj>(t))dt,  W  Il>  E  y  (25) 

J  —  OO 

with  y  =  {V>  e  C0(— oo,  +00;  V)  and  rj>’  E  L3(— 00,  +00;  H)}. 

Here  Cq{— 00,  +00;  V)  denotes  the  set  of  continuous  V-valued  functions  that  vanish  outside 
a  compact  time  interval. 
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In  [10,  9],  it  has  been  proven  that  a  necessary  and  sufficient  condition  for  a  T-periodic 
solution  to  be  a  weak  solution  in  (—00, -foo)  is  that  it  must  be  a  weak  solution  in  the 
interval  [0,T].  The  regularity  of  periodic  solutions  can  be  obtained  by  noting  that  3fx  6 
(— oo,  -foo)  such  that  ue(x ,  t\)  £  V.  We  use  this  as  the  initial  value  and  study  the  regularity 
in  [<i,  oo)  by  standard  methods  [15].  Repeating  this  in  each  interval  we  deduce  the  regularity 
of  the  periodic  solution  in  the  entire  interval  (— oo,+oo)  as  ue(t)  £  Z/°°(  — oo, -f  oo;  V)  fl 
L2loc{  — °Oj  +°o;  D(A)).  Let  us  now  state  the  main  theorem: 

Theorem  5.1  Let  f(t)  6  L2loc{  —  oo,  -foo;  H)  be  given  and  ue(t)  £  L°°{  —  oo,  +oo;  V)  fl 
LU  — oo,  -foo;  Z)(>1))  be  the  T-time  periodic  solution  of  the  regularized  Navier-Stokes  equa¬ 
tions.  Then  as  e  — »  0,  ut(t)  converges  to  a  limit  u  £  L°°(— oo,  -foo;  V')TlLfoc(— oo,  -foo;  D(A\)) 
which  is  a  T -time  periodic  solution  of  the  Navier-Stokes  equations 

Proof:  Since  ut(t)  is  a  T-periodic  solution  which  satisfies  (25),  for  any  k  £  Z  we  can  choose 
a  test  function  v(f)  =  V’oCO  so  that  it  vanishes  outside  the  compact  interval  (kT,(k- f  1)T). 
Thus, 

/■(Jfc-t-l)T 

{-(“<(*)»  V>o(0)  +  e(Aue(f),  AV>0(*))  +  i/(Vue(t),  Vv>0(t)) 

J  kT 

r(k+l)T 

-ffe(ue(t),«e(t),V>0(t))}dt  =  /  Vk£Z, 

JkT 

with  ut(t)  £  2z°°(fcT,  ( k  -f  1)T;  V )  fl  L2(kT,  ( k  -f  1)T;  D(A)).  In  particular  if  we  let  k  =  0,  we 
obtain  the  strong  solution  in  the  sub-interval  (0,T). 

Let  {uefc}  be  a  sequence  of  such  solutions  with  e*  =  1/k.  Then, 

/  {-(«€„(*),  V’o(O)  +  £fc  AV’oC*))  -f  i/(VutM(t),Vi/>0(t)) 

Jo 

+  K'iek(t)>U'k(t)>'f>o(t))}dt=  I  (/(*)>  O(0)^i 

Jo 

v  V’neCCo.Tjy)  ^0eL2(o,T-,H). 

From  Lemma  3.2,  we  have  the  strong  solution  utfc  £  T2(0,  T;  D(Ai))  D  L°°(0, T;  V")  with 
bounds  uniform  in  e*.  Thus,  by  applying  the  convergence  Proposition  2.1  we  obtain 

JO 
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+b(u(t),u(t),ij>0(t))}dt=  f  (f(t),ij>0(t))dt, 

Jo 

V  4oeC{0,T-,V)  i>'0el2(v,t-h). 

Here 

£k  f  (Aueji,Ai/> 0)dt — »0 
Jo 

as  e*  — ►  0. 

Now,  by  periodicity  of  the  strong  solution  we  get 

r{k+l)T 

{-(.(«),  v-i(t))  +  KVuti),  V*4(t)) 

JkT 

r(k+l)T 

+  Ku(t),u(t),r/>k(t))}dt  =  (f(t),1>k(t))dt,  (26) 

for  each  sub-interval  ( kT,{k  +  1)T).  Here  we  denote  V’t(i)  =  V>o(*  —  &T). 

The  above  result  implies  that  u(f)  €  L°°(kT ,  (A:  +  1)T;  y)  fl  L2{kT ,  (A:  -f  1)T;  D(Ai))  is  a 
strong  solution  in  each  sub-interval.  Hence,  by  choosing  a  suitable  i/>  (t)  and  summing  over 
A:,  we  get  (25)  with  e  =  0.  This  completes  the  proof. 

□ 


6  Upper  Semicontinuity  of  Attractors 

In  this  section  we  will  show  that  the  compact  global  attractor  Ae  associated  with  the 
semigroup  Wt(t,  •)  converges  to  the  compact  global  attractor  A  associated  with  the  conven¬ 
tional  Navier-Stokes  system.  We  prove  this  result  by  establishing  the  upper  semicontinuity 
of  Ae  at  e  =  0.  In  this  section  f  E  H  and  is  independent  of  time. 

We  begin  by  showing  that  W€(t,  •)  admits  a  compact  attractor  Ae  in  H  for  each  e  >  0. 
According  to  Ladyzhenskaya  [6],  the  formal  definition  of  a  global  attractor  is  as  follows:  A 
set  At  C  H  is  called  the  global  attractor  for  the  semigroup  Wt(t ,  •)  if 

(i)  A«  is  bounded  in  H\ 

(ii)  A€  is  an  invariant  set,  Wt(t,  0;A«)  =  A„  Vt  €  JR; 

(iii)  At  is  a  compact  set  that  attracts  the  bounded  sets  of  H. 
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For  a  study  on  the  global  attractor  to  the  Navier-Stokes  equations  the  reader  is  refered  to 
Constantin,  Foias  and  Temam  [3]  and  Temam  [16]. 

It  follows  from  Corollary  3.1  that  for  /  £  H  we  have 

l«.(f)|J  <  M0)|  V-  +  pi  [  1  -  e-“'  ],  t  >  0,  (27) 

where  a  =  vp,\  and  p0  =  |/|/a.  Note  that  both  a  and  p0  are  independent  of  e.  Hence,  for 
any  ball  #Ro  =  {«<(°)  £  l«*(°)l  <  Bo},  there  is  a  ball  BTq  in  H  centered  at  origin  with 

radius  r0  >  p0  ( R0  >  r0)  such  that 

W.(t,0;B„')cB„,  fort  >(„(£«,)=!  In  ^4- 

a  ro  -  Po 

The  ball  Bro  is  said  to  be  exponentially  absorbing  and  invariant  under  the  action  of  the 
map  We(i,0;). 

Recalling  the  inequality  (14),  we  have 

||«,|’  +  2e||«,f  +  <  \ll.  (28) 

Integrating  from  t  to  t  +  1,  we  obtain  for  tt0  G  Brq 

J‘+'  |V«,| \*ir  <  i  (rj  +  V  t  >  to (B^).  (29) 

Now,  we  recall  the  inequality  (17), 

iv«e|2  +  i/Mi|Vue|3<!^. 

CLT  V 

We  drop  the  positive  terms  vp.\ |Vu€j3  and  integrate  with  respect  to  r  from  s  to  t  -f  1  to 
obtain 

|v«.(t  +  l)p  <  |Vu.(j)|J  +  l^(t  +  1  -  *). 

We  then  integrate  the  above  inequality  with  respect  to  s  from  t  to  t  +  1  to  get 

|V«.(1  +  1)P  <  jf'+I  |v«.(5)|J<fa  + 
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It  follows  from  (29)  that 

|Vue(t)|2  <  l-  ^  ^  =Rl  V  t  >  t0( Brq )  +  1. 

This  means  for  €  B Ro  and  any  r  >  tQ(BR<>)  +  1,  we  have  ue(r)  €  Br v  =  {«t(r)  G 
V ;  |Vtit(r)|  <  i?x}.  That  is, 

We(r,0;  B^)  C  Brv,  for  t  >  t0 (B^)  +  1.  (30) 

From  Corollary  3.2,  it  is  easy  to  show 

|V«e(t)|3  <  |Vut(r)|2e-“(‘-r)  +  p\  (l  -  e~a^)  ,  Vi  >  t  >  to(B*)  +  1, 

where  a  =  j/pi  and  p\  =  |/|a/t/a.  Again,  both  a  and  px  are  independent  of  e.  Hence,  for 
any  ball  BRl,  there  exists  a  ball  Bn  in  V"  centered  at  origin  with  radius  rx  >  pi  ( Ri  >  rx) 
such  that 

We(t, r;  B* )  C  Bn ,  for  i  >  f^B*,)  =  t^B*)  +  l  +  -ln  ^4-  (31) 

a  rx  —  Pi 

The  ball  Bri  is  said  to  be  exponentially  absorbing  and  invariant  under  the  action  of  the 
map  Wt(t,T\-).  Combining  results  in  (30)  and  (31)  and  applying  the  semigroup  property 
gives  us 

=  W«(f,r;  W^t.OjBho)),  T>t0(BRo)  +  l 
=  W€(t,T-,BRl)cBri,  Vt>t1(BHo). 

That  is,  ball  Bn  of  radius  rx  in  V  will  absorb  any  bounded  set  Br 0  in  H.  Since  the 
embedding  of  V  in  H  is  compact,  we  deduce  that  ^*(^0;  •)  maps  a  bounded  set  in  H  into 
a  compact  set  in  H.  In  addition,  the  operators  Wt(t,  0;  •)  are  uniformly  compact  for 
t  >  fx( Br 0).  That  is, 

(J  Wt(t,Q\  BRo) 

»>‘i 

is  relatively  compact  in  H.  Due  to  a  theorem  from  Temam  [16,  Theorem  1.1],  there  exists  a 
compact  attractor  A*  which  attracts  every  bounded  sets  in  H.  A,  is  the  global  attractor  for 
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operators  Wt(t,  0;  •)  and  it  is  also  the  oi-limit  set  of  the  absorbing  set  Bri  (i.e.  Ae  =  uj(BTi)). 
This  means  if  we  denote  W€{t,  0;  BTl)  =  BTl(t)  then 

A.=  f|  c‘[UBrM 

T> 0  \t>r 

Note  that  the  global  attractor  Ae  must  be  contained  in  the  absorbing  balls  in  H  and  V: 

Ae  C  BronBri. 

Notice  that  all  the  above  bounds  are  independent  of  e.  In  particular,  for  e  =  0,  we  obtain 
existence  of  the  compact  attractor  A  for  the  conventional  Navier-Stokes  equations.  We  now 
prove  the  following  theorem: 

Theorem  6.1  For  e  >0,  We  admits  a  compact  attractor  Ae  which  attracts  bounded  sets  of 
H  and  is  contained  in  the  absorbing  balls  Bro  fl  Bri ,  where  ro  and  ri  are  independent  of  e. 
Moreover,  dff(Ae,  A)  — ♦  0  as  e  — *  0. 

Proof:  We  need  the  following  result  [4,  Lemma  2.1]: 

Let  X  be  a  Banach  space  and  T(t),  t  >  0,  a  semigroup  on  X.  Let  Thft )  be  an  approximate 
semigroup  (depending  on  a  parameter  h  >  0)  to  the  semigroup  T(t).  For  S  >  0,  let  8) 
denote  the  ^-neighborhood  of  a  bounded  set  B  €  X  which  is  the  union  of  open  balls  of  radius 
8  centered  on  B. 

Proposition  6.1  Let  B  €  X  be  a  bounded  set  and  Uq,U\  be  two  open  sets  such  that 
J\f(B,d0)  C  UQ,  Af(B,di)  C  U\  for  some  d0,di  >  0.  If 

(i)  B  attracts  Uo  under  T{t )  and 

(ii)  Thft)  approximates  T(t)  on  U\  uniformly  on  compact  sets  of  [to,  oo),  with  t0  >  0, 
then  for  any  8  >  0,  there  are  h0  >  0  and  r0  >  such  that,  for  0  <  h  <  h0,  for  t  >  t0, 

Th(t)(U0nUi)cM(B,6). 


□ 
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We  will  now  prove  the  semicontinuity  property  for  the  global  attractor  Ae.  This  can 
be  proven  by  showing  the  hypotheses  of  above  proposition  are  satisfied  by  !F(t,  (),•)  and 
We(t,0;-)  for  e  small  enough.  Clearly,  it  is  sufficient  to  show  that  the  ^-neighborhood  of 
attractor  A  is  an  absorbing  set  and  that  Wt(t,  0;  •)  approximates  W(t,  0;  •)  on  Br2  =  (u£(0)  G 
V]  |Vue(0)|  <  R2}  uniformly  on  compact  sets  of  [0,oo).  Let  us  prove  this  in  two  steps: 
First  step: 

Let  Af(A,8)  be  the  8- neighborhood  of  the  attractor  A.  Since  A  is  a  global  attractor,  for 
any  bounded  set  Brq  =  {“(0)  £  H\  |u(0)|  <  -Ro}  C  H 

dH(W{t,  0;  Bib),  A)  — »  0  as  t  ->  00. 

Here  dH(A ,  B )  is  the  semidistance  of  two  subsets  A,  B  of  H : 

dn(A,B)  =  sup  inf  d(x,y). 
x  eA  v£B 

Thus,  there  exists  8  >  0  and  t>ts  such  that 

for  t  >  ts- 


This  implies 


W{t,0\BRo)  for  t  >  is- 


This  shows  that  ^(A,^)  is  an  absorbing  set. 

Second  step: 

We  want  to  show  Wt(t, 0;  •)  approximates  W{t,  0;  •)  on  BRa  uniformly  on  compact  sets  of 
[0, 00).  By  subtracting  (13)  from  (12),  we  obtain  for  w  =  ue  -  u 
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Since  D(A1!2)  =  V  =  D(Ay),  the  second  term  of  left  hand  side  of  equation  (32)  can  be 
written  as 

e{Aut,w)  =  e(Aut,  ue)  —  e(Aue,  u) 

=  e\A1'2ut\2-e(A1'2ue,A1'2u) 

=  e\Ax,2ut\2  -  e(A1/2ue>  ^4iu). 

This  leads  to 

-Tr|w|2  +  e|>l1/2«£|2  +  i/|Vto|2  =  e{Ax^2ue,Axu)  +  b(tv,w,ve). 

z  at 

By  applying  Young’s  inequality,  we  have 

<£|^,/!»,|  |A,»| 

The  trilinear  term  can  be  estimated  as 

|6(io,tu,ue)|  <  Ci  (to  I1/2 1  Vto  p^ltD  |1|/,2|Vto  ^^IVujI 

<  Ci|w|  jVw  |  | V ucJ 

<^iv«i2+^r,v««|a. 

By  combining  all  of  the  above  estimates  in  (32),  we  get 

4 1®  I2  +  e \Ail2ut\2  +  v\ Vw  |2  <  6 1 Axu  I2  +  ^|to  |2| Vw  |2. 

at  v 

We  can  drop  the  positive  terms  e\Ax^2vit\2  and  t/|Vu>|2  to  obtain  the  following  differential 
inequality 

£W’<€|A1«|1  +  d|„|’|v,,|’.  (33) 

Using  Corollary  3.2  we  get, 

|V»,(()|J  <  |V»«(0)|V«‘  +  />!(!-  «-“')• 
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This  implies  that  any  ball  Br7  —  {ue(0)  6  V]  |Vue(0)|  <  R2}  in.  V  with  radius 
satisfy 


W£(t,0]BR3)  C  Br3,  for  t  >  0. 


This  means  if  ue(0)  (E  Br7,  then  We(t,  0;  ue(0))  is  defined  and  belongs  to  Br7  for 
ball  Br7  is  therefore  invariant  under  the  map  We(t,0;  •).  That  is  ,  we  have 


|Vue(0|  <  R2,  t>  0,  for  »e(0)  €  Br7. 


Let  uc(0)  6  Br7  then  it  follows  from  (33)  that 

d 


^M2<e|A«|2  +  ^|«,|2,  <>0. 


From  the  standard  Gronwall  Lemma,  (35)  gives 

M*)|2  <  |w(0)|2exp  (j  —————dT ^ 

l^iu(s)|2exp  ( 

rs) 

=  e  /  |^i«(s)|2  exp 
Jo 

r«?«(*-) 

* 

ds 

ds, 


where  t»(0)  =  ue(0)  -  u(0)  =  0. 

Let  us  now  consider  t  €  [0,T]  such  that  0  <  3  <  t  <  T.  This  gives 

\w(t)\2  <  e  exp  |v4xu(5)|2ds. 

By  applying  (ii)  of  Corollary  3.2  (with  e  =  0),  we  obtain 

M0I2  <  C-  exp  {^-Tj  (r22  +  ,  t  €  [0,  T). 

This  means  for  ue(0)  E  Br7,  then 


with 


|We(t,0;ue(0))-^(t,0; 


^m\2H<9(e,i'tR2,f,T), 


s(e,  v,  R3,f,T)  =  e 


exp 


'■2  >  Pi  will 

t  >  0.  The 

(34) 

(35) 
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and 

lim  g  (e,  u,R3,f,T)  =  0. 

Hence  We(t)  approximates  W(t)  on  Br3  uniformly  on  [0,  T].  According  to  the  Propo¬ 
sition  6.i,  then  for  any  6  >  0,  there  are  eo  >  0  and  To  >  0  such  that 

We(t)  ( Br,  n  BRl )  C  Af(A,  6 ),  for  0  <  e  <  c0,  t  >  r0. 

Since  the  attractor  Ae  is  contained  in  Brc  fl  Br3  ,  we  have 

W,(t)(Ae)  C  Af(A,  8),  for  0  <  e  <  e0,  t  >  r0. 

Since  Ae  is  an  invariant  set,  we  deduce  that 

Ae  C  Af(A,8),  for  0  <  e  <  e0,  t  >  r0. 

Since  8  is  arbitrary,  we  obtain  the  upper  semicontinuity  of  Ae  at  e  =  0: 

dji{ Ae,  A)  — ♦  0,  as  e  — >  0. 

□ 
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